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We probe the pure Einstein's gravity contributions to the second-order density power spectrum. 
In the small-scale, we discover that the Einstein's gravity contribution is negligibly small. This 
guarantees that Newton's gravity is sufficient to handle the baryon acoustic oscillation scale. In the 
large scale, however, we discover that the Einstein's gravity contribution to the second-order power 
spectrum dominates the linear-order power spectrum. Thus, pure Einstein gravity contribution 
appearing in the third-order perturbation leads to an infrared divergence in the power spectrum. 
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The large-scale cosmological density power spectrum 
can be regarded as one of the main pillars of modern 
cosmology where theories meet with observations. Re- 
cent discovery of the baryon acoustic oscillation (BAO) 
near 100/i _1 Mpc scale has spurred renewed interests 
in cosmology community on the importance of detailed 
theoretical studies of the physics behind the BAO in- 
cluding the initial spectrum and the nonlinear process; 
h is a Hubble constant in the unit of 100km/ sec/Mpc. 
The BAO provides an important distance scale of sound 
horizon size at the baryon decoupling epoch, which is 
completely set by the anisotropic angular power spec- 
trum of cosmic microwave background radiation (CMB), 
and hence, can be used as a standard ruler. This stan- 
dard ruler encoded in the power spectrum provides the 
measurement of the angular diameter distance and the 
Hubble expansion rate by which we can constrain the 
properties of dark energy. Recent studies show that 
nonlinear (next-to-leading order) processes in the power 
spectrum is important in theoretical analysis of the BAO 
phenomenon especially in the context of planned precise 
observations of the high redshift galaxy surveys [2| • 

Up until now, however, all studies of the cosmologi- 
cal power spectra are treated in the context of Newton's 
gravity. The linear-order cosmological perturbation was 
first handled in Einstein's gravity by Lifshitz in 1946 Q. 
Later a Newtonian study was made by Bonnor in 1957 
|4| , and it was shown that Newtonian density perturba- 
tion equation coincides exactly with the relativistic result 
in the zero-pressure limit. We have recently shown in Q 
that the density perturbation equation in Einstein's grav- 
ity coincides exactly with the previously known result in 
Newton's gravity Q even to the second-order: we call it 
a relativistic/Newtonian correspondence to the second- 
order perturbation. This is a nontrivial result which is 
available in the temporal comoving gauge with suitable 



identifications of relativistic metric and matter variables. 
In the zero-pressure case the Newtonian hydrodynamics 
is known to be closed to the second order in nonlinearity 
|6j. Thus, any nonvanishing third-order terms in Ein- 
stein's gravity can be naturally identified as pure general 
relativistic contributions. Again with the same comoving 
gauge and a suitable choice of variables recently we have 
derived the densityperturbation equation with the third- 
order corrections [7J ; all variables in our gauge condition 
are naturally gauge- invariant. 

The pure Einstein's gravity third-order terms con- 
tribute to the next-to-leading-order (second-order) den- 
sity power spectrum. In this Letter, we present con- 
tribution of the pure Einstein's gravity corrections to 
density power spectrum. Our results reveal cosmologi- 
cally important discoveries in both the small scale and 
the large-scale. Particularly surprising is the unexpected 
surfacing of an infrared divergence due to the general 
relativistic nonlinear effect in the large scale which po- 
tentially demands a new revised understanding of the 
current standard paradigm of the large-scale structure 
formation mainly based on the linear perturbation the- 
ory. 

We expand the density fluctuation as S = S± + 62 + S3 + 
■ ■ ■ where S(x,i) = Sg(x,t)/g(t) is the relative density 
fluctuation. The density power spectrum becomes 
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where A: is a wave number in Fourier space and () in- 
dicates phase space averaging; (Re(Sld2)) vanishes as- 
suming the random phase Q; a superscript * indicates 
a complex conjugation. The density power spectrum is 
introduced as (<5(k)<5(k')) = (2n) 3 S D (k + k')P(k) with 
k = |k| = |k'| and 80 a Dirac delta-function, [9j. The 
P\i is the linear-order power spectrum, and P22 + -Pi 3 is 
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the second-order (next-to-leading order) power spectrum. 
We notice that the third-order pure general relativistic 
correction terms contribute to P13; thus we decompose it 
as P13 = Pi 3 ,Nowton + -Pi3,Einstom ■ The density and veloc- 



ity power spectra valid to the next-to-leading-order are 
presented in [Io| . An i nteg ration of the density power 
spectrum in Eq. (19) of [l(| gives 
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where r = k'/k and x = (k • \t')/(kk'); l/i H = d/(kc) 
is the ratio between a scale £ = a/k and the horizon 
scale £h = c/H with 7? = a/a, and a(i) the cosmic scale 
factor. The second-order power spectrum in Eq. @ is de- 
rived in Einstein's gravity; as we have shown that results 
in Einstein's gravity coincide exactly with the Newton's 
one to the second-order perturbations, the power spec- 
trum in Eq. © without Pi3,Einstoin is exactly valid in 
Newton's gravity. The P22 and P13, Newton were presented 
in [ll| in the Newtonian context. The pure Einstein's 
gravity contribution P^Einstcin is the new result in this 
work; in Eq. ^ we assumes flat matter dominated back- 
ground. This pure Einstein's gravity contribution to the 
one-loop corrected power spectrum is hitherto unknown 
in the literature based on Newton's gravity [?|. We note 
that the pure Einstein's gravity contribution is multiplied 
by a (£/£h) 2 factor, thus suppressed far inside the hori- 
zon. The relativistic/Newtonian correspondence to the 
second order is valid in the zero-pressure situation with- 
out rotation, no gravitational waves, and spatially flat 
background, but valid in the presence of the cosmologi- 
cal constant; for the general cases, see 

In Figure [1] we present integration of equation by 
using a realistic linear power spectrum from a concor- 
dence cosmology of ACDM Universe. We calculate the 
linear power spectrum using the CAMB [13J code with 
the maximum likelihood cosmological parameters given 
in Table 1 of [HJ ("WMAP+BAO+SN"). Since equa- 
tion |2]) is valid for the spatially flat, matter dominated 
universe, we calculate the second (next-to-leading) or- 
der power spectrum at higher redshift (z — 6) when 
our Universe is well approximated by flat, matter dom- 
inated universe. When calculating the integration, we 
set the integration lower bound as r m - m — 10 -2 fcy/fc, 
where kn = 1/£h is the wave number corresponding to 
the horizon scale, and upper bound as r max = 10 3 /fc. 



Note that, unlike the Newtonian correction terms, P22 
and P13, Newton, P13, Einstein is sensitive to r max , and loga- 
rithmically diverges. It is a problem that the amplitude 
of the Einstein term contains a logarithmically divergent 
integral; a logarithmic divergence is also encountered in 
[l5l ] in considering the backreaction of inhomogeneity on 
the cosmological expansion [16]. 

In the small-scale it was found in 0] that the lead- 
ing order of P22 cancels exactly with the leading order 
of P13, Newton- Our result in Figure [1] shows that despite 
such a cancelation of the leading-order Newtonian contri- 
butions, the next leading order Newtonian contribution 
is still bigger than the pure general relativistic contri- 
bution P13, Einstein- We find that Pi3,Emstem is roughly 
1% of P22 + P13, Newton at fc = 0.1 h/Mpc, and it be- 
comes smaller and smaller as k increases. This may be a 
good news to the cosmology community based on New- 
ton's gravity because our result guarantees to use New- 
ton's gravity in handling the weakly nonlinear processes 
far inside the horizon including the BAO scale. In the 
large scale, however, Figure [1] reveals a completely un- 
expected surprising result: the pure general relativistic 
contribution to the second-order power spectrum domi- 
nates over the linear-order relativistic/Newtonian power 
spectrum. We emphasize that previous nonlinear per- 
turbation studies based on Newton's gravity have been 
mainly concerned with the small-scale effect where indeed 
the fluctuations grow from linear to nonlinear phases. We 
will discuss implications of this new discovery in the fol- 
lowing. 

Figure Q] apparently shows that Pi3,Emstein is bigger 
than even the linear-order density power spectrum Pn 
in the large scale; Pi3,Einstein is negative in general. It is 
not necessarily a problem that the Einstein term dom- 
inates at the smallest wavenumber k; this also happens 
for instance when a fluctuation distribution with no (very 
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FIG. 1: Second-order power spectrum and the contribution from each component of Eq. (2) at z = 6. Note that we take the 
absolute value for negative terms, and show with dashed lines. Vertical dotted line shows the wavenumber corresponding to the 
current comoving horizon (fcn). We use kmin = W~ 2 kH and k m ax = 10 3 to evaluate the integration in Eq. (2). [color online] 



little) large scale power at higher order grows a fc 4 tail at 
small k [la, E3- It i s > however, a problem that the domi- 
nant contribution is negative, which is a sign that pertur- 
bation theory is at least incomplete; the full power spec- 
trum is of course positive definite Jig]. At its face value 
our result implies that as the second-order effect is larger 
than the linear one the nonlinear effect of Einstein's grav- 
ity leads to breakdown of the perturbation theory in the 
large scale but far inside the horizon scale. The situa- 
tion is particularly awkward because in the equation level 
the third-order pure Einstein's gravity correction terms 
are smaller than the second-order relativistic/Newtonian 
terms by a factor (p v ~ <5$/c 2 , see Eqs. (l)-(3) in [lol ]; 
<5$ is the Newtonian gravitational potential which is 
generally quite small, 6$/c 2 ~ 10~ 7 - 10~ 5 . Thus, 
we naturally anticipate the perturbatively derived third- 
order solutions from the relativistic/Newtonian second- 
order equations [these are smaller than the second-order 
terms by a factor 5 ~ (£h/£) 2 8§] are bigger than the 
pure Einstein's gravity contributions inside the horizon, 
see Eq. ([2]). That is, in the large-scale limit (A; — > 
and r — > oo) a naive examination of Eq. @ shows 
that both Pi3,Newton and Pi3,Einstein have fc-dependence 
proportional to k~ 2 \5t(k, t)\ 2 J °° dk'(. . . ). However, an 
asymptotic expansion of Pi3,Newton shows that both the 
leading-order term (— 42r ) and even the next-to-leading- 
order term (100r 2 ) cancel away, whereas no such cance- 
lation occurs in Pi3,Einstcin- This explains the diverging 
fc-dependence of Pi3,Einstcin m Figure [T] compared with 



the behavior of Pi3,Newton- The divergence of Pi3,Einstein 
is not a real problem as the density fluctuation is de- 
scribed by V = (k 3 /2n 2 )P; as we have As, Einstein oc fe -1 , 
~Pi3, Einstein is convergent in that limit [lq |. 

A similar naive examination of Eq. Q in the small 
scale leads us to expect that P13, Einstein could be com- 
parable to the Newtonian second-order contributions. 
That is, in the small-scale limit (k — > 00 and r — » 0) 
a naive examination of Eq. ^ shows that P^Newton has 
a fc-dependence proportional to fc 4 \Si (fc, t)| 2 J °° dk'(. . . ), 
whereas Pi3,Einstein has a fc-dependence proportional to 
fc 2 \Si(k, t)\ 2 dk'(. . . ). However, as the leading-order 
terms in Pi3,Newton cancel exactly with the one in P22 [H, 
we have P22 + P13, Newton has the same fc-dependence as 
the leading-order term (— 45/(4r 2 )) in P^Einstein- Thus, 
we naturally anticipate that P13, Einstein could be compa- 
rable with the relativistic/Newtonian contribution to the 
second-order power spectrum in the small scale. How- 
ever, an asymptotic expansion of Pi3,Einstein shows that 
both the leading-order term (— 45/(4r 2 )) and even the 
next-to-leading-order term (—21) cancel away, whereas 
no such cancelation occurs in Pi3,Newton- This explains 
why Pi3,Einstein is far smaller than the Newtonian contri- 
butions in the small scale. 

As a matter of fact, we can hardly be able to accept 
the result that the Einstein's gravity leads to breakdown 
of the perturbation expansion in the scale where per- 
turbation amplitude is known to be near linear. Being 
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confronted by this unexpected situation we have gone 
through all the algebra several times without finding a 
mistake in the analytic calculations leading to equation 
@; for convenience in the close examination of the al- 
gebra by the readers we made the whole calculation in a 
pdf file available in the web [HI]; the numerical code is 
also uploaded in the web . Besides any potential error 
made in our part, we point out that the infrared diver- 
gence could be due to an improper choice of our gauge 
condition. Our results are based on the temporal comov- 
ing gauge, setting T° = 0. The relativistic/Newtonian 
correspondence to the second order is available only un- 
der this gauge. We have continued to use the same gauge 
condition and identification of perturbation variables [7j ■ 
Whether we could find a nondivergent power spectrum 
in other gauge condition is an interesting possibility to 
be investigated in the future work. 

Recently, Losic and Unruh presented an intriguing pos- 
sibility of divergent behaviors of second-order perturba- 
tions in the quantum generation stage durin g in flation 
which might have a close relation to our result [2fJ . They 
argued that "a certain nonlocal measure of second-order 
metric and matter perturbations generically dominates in 
its amplitude compared to that of the linear-order pertur- 
bations during slow-roll inflation" . They further argued 
that "during slow roll, second-order fluctuations grew 
large for a class of inflationary models" and "nonlinear, 
and probably nonperturbative, gravitational effects dom- 
inate near slow-roll spacetime, and therefore the linear 
perturbation theory likely fails in those situations" . The 
linear power spectrum used in Figure [TJ is based on a near 
Harrison-Zel'dovich spectrum (21| which naturally arises 
from quantum fluctuations during the slow-roll inflation 
era (22j; this is often regarded as the major triumph of 
the inflation scenario in the early universe and helped to 
make inflation scenario a firm theoretical feature in the 



early universe despite its energy scale far beyond exper- 
imentally reachable range. Our result in this Letter is 
analogous to and consistent with the argument made by 
Losic and Unruh. Based on the near Harrison-Zel'dovich 
spectrum, we find a possible divergence in the large scale 
due to pure Einstein's gravity effect on the second-order 
power spectrum. One way out of this conundrum is to 
have a correct initial power spectrum generated in the 
inflation era including the role of the second or higher 
order perturbation theory during slow-roll inflation era. 

Divergent results found in the nonlinear perturbations 
in the large scale of density power spectrum here and 
in the seed generation stage of early inflation era found 
by Losic and Unruh, if confirmed to be correct, seriously 
challenge the currently accepted standard paradigm of 
theoretical cosmology concerning the cosmological struc- 
ture formation theory and physics in the early universe. 
A systematic investigation of the second-order perturba- 
tion theory during the quantum generation stage based 
on inflation is needed to resolve the issue raised in [2C| . 
Our result also suggests systematic investigations to be 
made concerning the role of third-order perturbations on 
the CMB temperature anisotropy power spectrum in the 
large scale. 
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